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A model of phase transitions with coupling between the order parameter and its gradient is 
proposed. This nonlinear model is shown to be suitable for the description of phase transitions 
accompanied by the formation of spatially inhomogeneous distributions of the order parameter. 
Exact solutions of the proposed model are obtained for the special cases which can be related to 
the spinodal decomposition or cosmological scenario. Such solutions manifest crucial dependence 
on the coupling constant. 
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Within the context of the general theory of phase tran- 
sitions, a system treated as a continuous medium is as- 
sumed to have a ground state which can always be de- 
scribed in terms of the order parameter. Usually, the 
ground state is described by a definite value of the or- 
der parameter that is independent of the point of the 
elastic continuum. The order parameter associated with 
the phase transitions can have various geometrical pre- 
sentations, for example, a scalar field in the case of con- 
densed matter [1 , or a fundamental scalar field in the 
field theory [2]. To introduce the order parameter that 
determines a stable state of the condensed matter, we 
have to consider probable deformations of the field dis- 
tribution, in particular, the disordered configuration of 
the ground state. Such distributions of the order param- 
eter can be observed experimentally, this issue requires 
an appropriate theoretical description. Today there ex- 
ist many phenomenological models which are declared as 
those that can give fairly general description of phase 
transitions and the behavior of spatial distributions of 
order parameters both before and after the phase transi- 
tion. To these well-known phenomenological models one 
can attribute the Landau theory of phase transitions pQ 
and the gradient theory of phase transitions [3], [4]. In 
the case of the first model, the free energy of the system 
under consideration is described in terms of the order pa- 
rameter while in the second case the description is given 
in terms of the order parameter gradient. Both models 
provide an efficient description of many details of phase 
transitions, in particular, the post-transition behavior of 
spatial distributions of the order parameter. 

The purpose of the present paper is to show that the 
standard model of phase transitions can be unified with 
the gradient theory by introducing in the standard model 
the coupling between the order parameter and its gradi- 
ent. In this case we get rid of the requirement to describe 
only phase transitions with homogeneous distributions of 
the order parameter. Namely, problems of this kind arise 
in modern physics, including the spinodal decomposition 
[5] and the alternative cosmological model that takes into 
account the formation of inhomogeneous distributions in 



the fundamental scalar field [2]. As it was shown in [5], 
the coefficients before the quadratic terms of the order 
parameter and its gradient can change their signs at the 
instability thresholds and thus some decomposition sce- 
nario becomes possible. This can result in the transition 
from the disordered state to modulated, patterned, or 
ordered-patterned states. In particular, in the case of 
cosmological scenario, it would be highly desirable to de- 
scribe the generation of spatially inhomogeneous states 
associated with the specific behavior of the coefficients 
in terms of the expansion of the model potentials. When 
the system described by a scalar spatially-dependent or- 
der parameter is rapidly quenching from a homogeneous 
high-temperature phase to the phase-coexistence region, 
we have the spinodal decomposition. The decrease of 
the temperature can make the reason for the formation 
of a new bubble phase in the cosmological model. Usu- 
ally, the theoretical analysis of the spinodal decomposi- 
tion and spatial formation of the fundamental scalar field 
is based on the standard model. In this case, however, 
the presentation for the free energy contains no terms 
that would take into account the interaction between the 
order parameter and its gradient though this interaction 
can lead to the restriction of the governing parameter. 

Any deformation of the order parameter requires addi- 
tional energy and thus imposes restrictions on probable 
deformations. In other words, not all deformations can be 
realized in system. In the general case the requirement of 
consistency of the known solutions of the standard model 
with the generalized description can be provided by an 
appropriate choice of the coupling between the order pa- 
rameter and its gradient. As it was mentioned above, the 
phase transitions associated with the system with broken 
continuous symmetry can be described in terms of the 
relevant order parameter. In particular, according to the 
Landau theory, the free energy density can be presented 
in terms of the order parameter as given by 

F = C(V^(r)) 2 + W^(r)) (1) 

where (p(r) is the order parameter, W(cp(r)) describes 
the order-parameter dependence of the free energy that 
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is assumed to be known, and C is a constant. 

In the well-known standard model of phase transitions 
we have 

^(r)) = ~VW + JV(r). (2) 

With the dimensionless variable (/? 2 (r) = Bip ^ being in- 
troduced, the standard dimensionless free energy reduces 
to the form given by 

/ = Z 2 (V^(r)) 2 + ((<p(r)) 2 -l) 2 (3) 

or 

/ = ^ 2 (V^(r)) 2 + ((^(r)) 2 -lV 2 (r) (4) 

with the potential being written in the standard form 
and I 2 = ^ being the characteristic length. It should 
be noted that in this presentation the dimensionless free 
energy depends on the relation between A and B and 
its sign. Making use of this expression for the free en- 
ergy density in mean field approximation, we can find 
the spatial distribution of the order parameter and thus 
to describe the properties of the new states which can be 
formed after the phase transition. The first presentation 
of the free energy density generates the Euler-Lagrange 
equation for the stationary distribution of the order pa- 
rameter that is given by 

Z 2 (V)V(r)=2<p(r)(^(r)-l) (5) 

and has a topological soliton solution cp = tanh(j). In 
the second case we can obtain another form of the Euler- 
Lagrange equation, i.e., 

? 2 (V)V(r) = ^(r)(l-2^ 2 (r)) (6) 

The latter equation has a different solution, (p = 
sech( j), which contains no topological singularities. This 
example shows that different coefficients in the expansion 
of the free energy result in different spatial distributions 
of the order parameter which can be observed experimen- 
tally. 

In the case of a system with the gradient of the or- 
der parameter, it is likely to be reasonable to introduce, 
in the free-energy functional, a term responsible for the 
probable interaction between the order parameter and 
its gradient. This coupling can regularize possible per- 
turbations of the order parameter and thus confine the 
spatially inhomogeneous state of the system. This means 
that not any possible deformation of the order parame- 
ter can exist in the deformed matter. Any defect in- 
troduces additional energy and not any deformation can 
contribute to the spatial distribution of the order param- 
eter. 

The main idea of the present contribution is to gen- 
eralize the phenomenological theory of phase transitions 
by introducing the interaction between the order param- 
eter and its gradient. Such generalization is reasonable 



in view of both mathematical and physical arguments. 
With regard to the above assumptions, the free energy 
can be presented in the form 

/ = a(|^(r)|) 2 + 6((|^(r)|) 2 -l) 2 + 

W(r)(lMr)) 2 1 1 

where operator D = /V is introduced, with the length of 
change of the order parameter; Z, a, 6, and c are parame- 
ters describing the influence of the gradient of the order 
parameter, the order parameter, and the coupling be- 
tween the order parameter and the order parameter gra- 
dient. If the coefficient c is equal to zero, then we return 
to the standard theory of phase transitions or the the- 
ory of fundamental scalar field. Nevertheless, these con- 
siderations demonstrate the need to treat the variation 
problems for the functional in various representations. 
We prove that, as the relaxation parameter vanishes, the 
families of such fields with finite free energy are compact 
and have many solutions which describe spatially inho- 
mogeneous distributions of the order parameter and pos- 
sible topological structures of the new phase. Our proof 
is based on the analysis of the extremum of the free en- 
ergy density. The method also yields an Euler-Lagrange 
equation for the limit of the defect structure measure. 
This measure also satisfies the cancellation properties de- 
pending on its local regularity, which seems to indicate 
several levels of singularities in the limit. We show that 
the same solution can be obtained in other interpreta- 
tion. To do this, we have to analyze the behavior of the 
free energy density. The free energy minimum satisfies 
the Euler-Lagrange equation 

(a + ap 2 (r))D 2 tp(T) + c^r)(Zty(r)) 2 ^ 
+2Mr)(l-(^ 2 (r)) = [ ' 

Let us consider several possible solutions of the Euler- 
Lagrange equation for various combinations of coeffi- 
cients in the one-dimensional case, when Dcp = = 
(p. The Euler-Lagrange equation can be rewritten in a 
simple form, i.e., 

i3 + 5 (^f + %) = (9) 

where g(ip) = 2 an d h(ip) = 2b ^ c ~2 ~ If a new vari- 
able u{ip) = ((/3) 2 is introduced, then the equation thus 
obtained can be reduced to the Bernoulli equation, i.e., 

u(<p) + 2g(<p)u(<p) + h(<p) = (10) 

where u(ip) implies derivation with respect to cp. The 
general solution of the Bernoulli equation in our case has 
the form 

uM = ^T^ (s + 6(1 "^ )2) 

where s is the integration constant which in the general 
case can be put equal to zero. From this solution of the 
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Bernoulli equation we obtain that the spatial behavior of 
the order parameter is described by the relation 

which yields an expression for the spatial distribution of 
the order parameter given by 

x = (§)*Ar*M(§)'¥>) + (^Arsh((^ 

The behavior of solutions for different signs of the cou- 
pling constant is shown in Fig.l. Thus, we can obtain 
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FIG. 1: Different solutions for various values of the coupling 
constant 

from this solution all the well-known standard solutions. 
For c = 0, we find that 

x= 1 -( b -)iArth(^-) (14) 
la (p z — 1 

or tp = th(ffix that fully recovers the solution of the 
standard theory of phase transitions, or the behavior of 
the fundamental scalar field in the standard cosmolog- 
ical model. If the coupling constant becomes negative, 
then the solution crucially changes. In the case when the 
c = —a, we obtain a periodic solution <p = sin(^)^x. 
This solution is not trivial in the physical sense. Within 
the context of this solution we can propose a new sce- 
nario of phase transitions. Namely, we can predict the 



existence of various spatial distributions of the order pa- 
rameter after the formation of the new phase. In order 
to make such a prediction we have to know (or postulate) 
the temperature dependence of the coupling constant. If 
the coupling between the order parameter and its gradi- 
ent is included into the free energy with positive value, we 
can observe an inhomogeneous distribution of the order 
parameter or the fundamental scalar field without topo- 
logical singularities. If the coupling constant changes as 
temperature changes and if its temperature-dependence 
is critical, c = c(T — T s ), where T s is the new critical 
value of temperature for the coupling constant equal to 
zero, then we can observe solution of the standard model 
of phase transitions, tp = th(^)^x with the topological 
singularity. The further decrease of temperature can lead 
to the negative value of the coupling constant with the 
dependence c — c(T — T s ) = —a and the topological so- 
lution transforms to the periodic solution cp = sin(^)^x, 
i.e., a regular solution without any singularities. 

Assuming that the interaction between the order pa- 
rameter and its gradient induces selection of possible de- 
formations that depends on the coupling constant, we can 
obtain various spatial distributions of the order parame- 
ter. Such phase transitions can be observed in the exper- 
iments with decreasing temperature. This effect can be 
expected in the standard cosmological model, when we 
have evolution of the bubble formation (as a topological 
formation) to a new periodic phase of the fundamental 
scalar field without singularities. This effect should be 
observed in the experimental studies of the spinodal de- 
composition. 

To conclude, we note that the model of phase transi- 
tions with coupling between the order parameter and its 
gradient can describe a new scenario of the spatial distri- 
bution of the order parameter after the critical point of 
the phase transition. In this model we obtain an exact 
analytical solutions and predict the experimental obser- 
vation of the new behavior of the order parameter in the 
case of the spinodal decomposition or in the cosmological 
scenario in the field theory. It is shown that this nonlin- 
ear model can be used for the description of phase transi- 
tions accompanied by the formation of spatially inhomo- 
geneous distributions of the order parameter for various 
systems. Such solutions crucially depend on the coupling 
constant. In the case of its negative value we can expect 
that the distribution of the order parameter should be 
periodic, but in the case of a positive coupling constant, 
we obtain solutions with topological singularities. 
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